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ABSTRACT 


A  series  system  for  which  only  failed  components  are 
replaced  (or  repaired)  is  considered.  Nonf ailed 
components  are  in  a  "state-of-suspended  animation" 
while  the  system  is  down.  The  limiting  average  system 
up  time  is  computed  for  arbitrary  failure  and  repair 
distributions.  The  limiting  distribution  of  system 
up  time,  the  number  of  failures  of  component 
i  (i  =  1,2,  ...,  k)  and  the  down  time  of  component 
i  (i  *=  1,2,  ...,  k)  are  calculated.  The  asymptotic 
distribution  of  the  cost  of  repair  is  also  derived. 


AVAILABILITY  THEORY  FOR  MULTICOMPONENT  SYSTEMS 


by 

Richard  E.  Barlow  and  Frank  Proschan 

0.  INTRODUCTION  AND  SUMMARY 

In  this  paper  we  consider  one  of  the  most  basic  stochastic  models  in 
reliability  theory — namely  the  on-off  process  generated  by  failures  and  repairs 
of  components  in  a  series  system.  A  series  system  of  k  components  operates  if 
and  only  if  each  of  the  k  components  operates.  No  component  operates  while 
the  system  is  down.  Furthermore,  only  failed  components  are  repaired  and/or 
replaced;  repair  or  replacement  takes  a  random  time.  Repaired  components  are 
assumed  to  function  like  new  components. 

There  is  a  large  literature  dealing  with  availability,  the  probability  that 
the  system  is  functioning.  However  most  papers  assume  special  repair  or  failure 
distributions  (or  both).  An  example  of  this  type  is  the  paper  by  Gaver  (1963) 
and  the  paper  by  Obretenov,  Dimitrov  and  Uzunov  (1969).  Many  papers,  including 
those  just  named,  are  concerned  with  parallel  systems  with  independently  operating 
components  -  that  is  nonfailed  components  are  not  usually  in  a  state  of 
"suspended  animation"  during  repair  of  a  failed  component,  as  In  our  model.  The 
reader  may  consult  the  IEEE  Transactions  on  Reliability  as  well  as  the  Colloquium 
on  Reliability  Theory  held  at  Tihany,  Hungary,  16-19  September,  1969.  A  model 
more  general  than  ours  is  treated  by  M.  C.  Botez  (1969).  However  there  is  no 
overlap  with  our  results. 

Let  be  the  duration  of  the  rth  functioning  period  of  component  i  with 

distribution  F^  and  mean  y^  ,  i  «  1,2,  ...,  k  (i.e.,  time  to  failure  of  rth 
replacement  for  component  i  excluding  system  down  times) .  Let  D^r  be  the 
duration  of  the  repair  time  (or  down  time)  for  component  i  with  distribution 
G_^  and  mean  ,  i  =*  1,2,  . . . ,  k  .  We  assume  that  for  i  =  1 . k  , 


■~*f«  —  ,. 
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|X_,  l  and  |D,  |  are  mutually  independent  renewal  processes. 

r=l,2,...  *  ir  r=l,2, . . . 

A  typical  failure-repair  history  for  such  a  series  system  might  look  as 
follows: 


FIGURE  1:  HISTORY  OF  A  SERIES  SYSTEM 


In  Figure  1,  component  i  fails  at  time  X^  and  the  system  is  down 
hours.  The  system  again  operates  from  time  X^  +  to  time  X^  +  . 

Component  j  fails  at  time  and  is  replaced  by  time  X^  +  +  ®jl  * 

Let  C(t)  =  i  if  the  system  is  down  at  time  t  due  to  the  failure  of 
component  i  (i  =  1,2,  ...,  k)  .  Let  5(t)  =  0  if  the  system  is  operating  at 
time  t  .  We  are  interested  in  the  limiting  probability,  lim  P[£(t)  =  i]  , 

that  the  system  is  in  state  i  .  The  limiting  system  availability,  lim  P[£(t)  0]  , 

t-*®° 


is  of  special  interest. 

Since  only  failed  components  are  replaced  with  new  or  like-new  components, 
the  age  distribution  of  components  in  the  system  quickly  becomes  mathematically 
very  complicated.  The  process  (?(t)  ;  t  >_  0 }  has  in  fact  no  regeneration  points. 
It  is  remarkable,  however,  that  many  quantities  of  interest  are,  in  the  limit, 
mathematically  very  simple  and  depend  only  on  component  mean  lives  and  component 
mean  repair  times.  The  limiting  average  system  availability,  as  we  shall  prove  in 
Section  2,  is 
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(0.]) 


PU(u)  =  0]du  =  lj  1  +  l  vjyj  d=f  7t 0 


while 


(0.2) 


lim  |  f  PtC(u)  =  i]du  =  --  nQ  d=f  ...  ,  i  =  1,2 . k  . 

t-*”  J0  1 


These  formulae  are  true  for  arbitrary  failure  and  repair  distributions.  If 

lim  P[5(t)  =  i]  exists,  then  it  is  equal  to  v.  .  Although  (0.1)  is  well  known 
t-x»  1 

for  exponential  failure  and  exponential  repair  [cf.  U.S.  Naval  Weapons  Reliability 
Engineering  Handbook  (1968)],  a  rigorous  proof  seems  to  be  missing  for  the  general 
case.  A  heuristic  proof  for  (0.1)  was  offered  by  Bazovsky,  MacFarlane,  and 
Wunderman  (1962). 

Let  N^(t)  be  the  number  of  replacements  of  component  i  in  time  t  .  We 
show  that 


EN.(t)  tt0 


i  “  1,2,  ...»  k 


This  result  can  be  used  in  determining  spare  parts  requirements,  since  ir^t/p^ 

will  be,  approximately,  the  number  of  components  of  type  i  required  in  [ 0 , t ]  . 

In  Section  3  we  prove  the  asymptotic  normality  of  t  **[^(0  -  where 

-I 

mi  =  ^i^O  is’  approximately,  the  mean  time  between  failures  of  component  i  . 


We  also  show  that 


t  ^[N(t)  -  tm 


is  asymptotically  normal,  where 


N(t)  -  l  N  (t) 


4 


A  similar  result  hojr  for  U(t)  ,  system  up  time  in  [0,t] 


t  ^[U(t)  -  flgt] 

is  asymptotically  normaj  arui.s  are  used  to  calculate  repair  and  maintenance 


costs  in  Section  4. 
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The  following  theorem  and  its  corollary  can  be  proved  by  the  argument  in  Feller, 
Volume  I,  (1968),  p.  321. 

1.2  Theorem: 

(1.3)  (on)  **  l  (X  -  p)  -  N(0,t) 

r=l 


as  n  ->  ~  if  and  only  if 

(1.4)  (o2u  3n)  ’5[N(nt)  -  ni/y]  -►  N(0,t)  , 

where  [  J  means  greatest  integer  contained  within  the  brackets. 

1.3  Corollary: 

oo  2 

If  is  a  renewal  process  with  EXr  =  p  and  Var  X^.  =  a  ,  then 

r=l  1 

both  (1.3)  and  (1.4)  hold. 

Billingsley  (1968),  pp.  148-150,  and  Iglehart  and  Whitt  (1969)  generalize 
Theorem  1.2  to  Wiener  processes  on  [0,1]  . 


where  (N^(t)  >  t  _>  0}  is  the  renewal  counting 

00 

l*,ir}  '  However,  N.(*)  and  U(t)  are  not 
r=l  1 

U(t)  =  t  implies  N^[U(t)]  =  0  .  (We  assume  that  all  components  are  new  at 
t  -  0  for  definiteness.  However,  the  limiting  results  are  true  regardless  of 
the  initial  conditions.) 

2 . 1  Lenuna : 

If  0  <  ui  <  «  and  0  £  v  <  »  (i  =  1,2,  k)  ,  then 


process  associated  with 
independent,  since  in  particular 


(2.1) 


lim 

t-x» 


Ni[U(t,w) ,w ] 
U(t,w) 


3  .  S  •  j 


where  we  have  included  the  argument  w  e  fi  to  emphasize  that  all  random  variables 
are  defined  on  the  same  probability  space  (ft,A,P)  . 


Proof: 

Under  the  hypotheses  U(t,w)  »  almost  surely  as  t  ■+•  °°  .  Since  N^(t,w) 
and  U(t,w)  are  defined  on  the  same  probability  space,  (1.1)  implies  (2.1). || 


2.2  Theorem: 


If  0  <  pi  <  «  and  0  vj,  <  oo  (i  =  1,2,  . . . ,  k)  ,  then 


(2.2) 


def 


a.s. 


Proof : 


Note  that 


k  k  »,<« 

I  I  D.r<D(t)<  l  I  D 


i=l  r~-l 


i=l  r=l 


The  inequality  results  from  the  fact  that  the  system  may  be  down  at  time  t 
Since  U(t)  +  D(t)  =  t  , 


k  Ni(t)  N  JU(t)] 

1+  l  r-1—  l  D  ir~j7^ 
i=l  Ni(t)  r=l 


Ni(t) 


By  the  strong  law. 


—  Id.  +  v .  as  t  ■ 
/.s  ,  ir  a.s.  i 


N.(t)  r=l 
1 


By  Lemma  2.1, 


N.[U(t)J  x 

— — — -r —  ->  —  as  t  -*■  00  .  Hence, 

U(t)  a.s.  p 


,  ,  U(t)  1 

a.  — > — r  *  "o  • 

i+  I  - 

i=i  “i 


The  reverse  inequality  is  proved  similarly. 


2.3  Corollary: 


Under  the  same  conditions, 


EU  (t) 

— £ —  ->  7iQ  as  t  -> 


Proof : 


Since  -  -  £  1  and  ------  -»■  rn  ,  it  follows  by  the  Lebesgue  dominated 

t  ’  t  3  •  S  •  U 

convergence  theorem  that 


EU(t) 

— £ —  ^  as  t  00  ’ 
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2.4  Corollary: 


Let  D^(t)  be  t*ie  down  time  for  component  i  in  [0,t] 


the  same  hypotheses, 


(2.3) 


lim 


D5(t) 


t  a.s.  p  0 
t-*°  i 


Proof : 


Note  that 


N.(t)-1 


Ni(t) 


l  D.  <  D . (t)  <  I  D.  , 
lr  —  l  —  lr 

r-1  r=l 


so  that 


yt* :  i 


^(t) 


Nj.(t) 


l  D 


r-1 


lr 


NilU(t)]  ^ 


U(t) 


Hence, 


Di<0  \ 
llm  —  <  —  »„ 

,  t  =  p,  0 

t-x"  i 


by  the  strong  law,  Lemma  2.1  and  Theorem  2.2. 


The  reverse  inequality  is  proved  similarly. | 


Of  course,  it  also  follows  from  Corollary  2.4  that 


k  v. 


t-*» 


(JO  =  w  y  _i 
t  a.s.  0  p. 


i=l  H1 
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■  iv ss.gKas',  t. v,\ 


-  ~.  ^?pss5^>j?5t>t^.>  ;«5  (iff^n3K,¥W^Q:  ■ 


Average  Availability 

We  call  T  *  j  P[£;(t)  =  0]dt  the  average  availability  in  [0,T]  .  It  is  a 
0 

well  known  property  of  stochastic  processes  that 


T 

(2.6)  T-1 


an  easy  consequence  of  Fubini's  Theorem.  It  follows  from  Corollary  2.3  that 


JpitM  -oidt-SMa. 


(2.7) 


lira  T 
T-k» 


P[C(t)  =  0]dt  = 


TI 


0  * 


If  lim  P[£(t)  =  0]  exists,  then  it  can  easily  be  shown  using  (2.7)  that 

t-*» 


(2.8) 


lira  PU(t)  «  0]  »  ti 

t-KO 


The  limit  in  (2.8)  will  not  always  exist  under  the  hypotheses  of  Theorem  2.2. 
Sufficient  conditions  for  (2.8),  for  example,  are  nonlattice  and  F^ 

exponential  (j  ?  i)  . 

Similarly,  if  P[4(t)  =  i]  exists,  then 


(2.9) 


lim  PU(t)  “  i]  -  v  . 

t-X»  1 


System  Mean  Time  Between  Failures 

Each  time  the  system  is  repaired,  the  time  until  next  failure  will  of  course 

depend  on  the  repair  history  of  each  component.  However,  the  average  of  successive 

up  times  will  converge  to  a  limit,  say  p  .  Likewise  the  average  of  successive 

down  times  will  converge  to  a  limit,  say  v  .  To  calculate  these  quantities,  let 
k 

N (t)  =  l  N. (t)  be  the  number  of  system  failures  in  [0,t]  . 
i=l 


tttl 1 1  Miai.' ■- r ■Vi-ir- r~ ■'  - - '< ^ 
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2.7  Theorem; 

If  0  <  <  00  and  0  <_  <  »  (i  «  1,2,  k)  ,  then  the  limiting  average 

of  system  up  times  will  be  a.s. 

(2.10)  P  -  1  /  l  ~ 

/  i=l  1 i 

while  the  limiting  average  of  system  down  times  will  be  a.s. 

k 

(2.11)  V  =  y  l  V  /w  . 

i=l  1 

Proof : 

The  average  of  system  up  times  in  fO,t]  will  be  approximately  U(t)/N(t)  . 
(The  error  will  go  to  0  a.s.  as  t  •*■  »  ,  as  in  previous  proofs.)  Since  by 

N.(t) 

Theorem  2.2,  lim  U(t)/t  =  ,  and  by  Corollary  2.5  lim — - —  =  tt n/y  , 

a.s.  u  t  a.s.  u  i 

t-*»  t-*00 

it  follows  that 

,  V  £  0(0/i(„  ”o/(”“iU)- 

Hence, 

lim  U(t)/N(t)  =  1  /  7  —  . 

.  a.s.  /  ,*•,  y. 

t-*»  /  i«l  i 


The  average  of  system  down  times  in  [0,t]  will  be  approximately 


N.(t) 

l  Dir/N(t)  . 
r*=l 


By  Corollary  2.5  and  the  strong  law. 


k  N.(t) 

l  e - 

i=l  N (t)  N  (t) 


iyt) 

l 


r=l 


-y 

a.s . 


V  d?f 

i=lWi 


where  now  y  is  defined  by  (2.10)  and  v  by  (2.11).  From  Theorem  2.7,  we  see 
that  (2.13)  is  the  analogue  of  (2.12). 
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3.  ASYMPTOTIC  DISTRIBUTIONS 


k 

To  obtain  the  asymptotic  distribution  of  N. (t)  ,  N(t)  **  £  N. (t)  ,  and 

1  i*l  1 

2  2 

U(t)  ,  we  must  also  specify  the  variances  **  Var  and  *  Var 
(i  =  1,2,  k)  .  Let  m^  =  vi^iTq^  for  i  ■  1,2,  . ..,  k  .  As  we  shall  showt 

,  is  approximately  the  mean  time  between  failures  of  component  i  .  Let  * 
denote  nortned  random  variables.  Although  normed  random  variables  will  have 
asymptotic  mean  0  ,  the  asymptotic  variance  is  not  necessarily  1.  In  particular, 
let 


(3.1)  N*(t)  -  t^Ct)  -  tm^1]  . 

We  first  state  our  main  results  before  presenting  proofs. 


3.1  Theorem: 

2  2 

If  0<p4<“>,0£Vi<®,0<o^<»,0<^T^<°o  for  i  «*  1,2,  . . . ,  k  , 

then 


[N*(t),N*(t) 


> 


is  asymptotically  (t  -*•  <*>)  multivariate  normal  with  nean  vector  0  and  variance- 
covariance  matrix 


t  -  (Vy) 


wnere 


(3.2) 


-1  -if  r  r  2  2  -3  ,  2  -ll  2-2  2-2  ,  .. 

=  “i  mj  [O  LVsVs  +  Vs  J  '  Viyi  “  Vj  Yj  J  (i  * 


-  “TV 


^  w  -  v  *-  » -^;  -,,-~  v:r' - *  , 
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(3.3) 

and 

and 

(3.4) 

3.2  Corollary: 

Under  the  conditions  of  Theorem  3.1, 

-  tm^1] 

is  asymptotically  (t  -*•  “)  N(0,1)  . 

3.3  Theorem: 

k  -  -1  k  -1 

Let  N(t)  *  £  N  (t)  and  m  »  £  m  .  Under  the  conditions  of 

i°l  1  i=l  1 


2  -3  2 

V  ay  ss  tn  w 

vii  i  i  wi 


2  22.  c  z  —l  V  z  z  - 

w  -  o  c  +  11.  I  t  y  +  y  l  0  V  U 

i  i  i  1  j.i  J  3  1  j*i  j  j  j 


2  -1 


2  2-3 


:i  ■  1  +  l  V?  •  ’ 


1  + 


Theorem  3.1, 

t  ^[N(t)  -  tra 


is  asymptotically  (t  ■*  «)  N(0,1)  ,  where 

..  2  -2  r  f,  ,2  2-3,  2  -ll 

(3.5)  v  *  Ti0m  l  ^(m  -  vi)  j  . 

\ 

3.4  Theorem: 

Let  U(t)  be  the  system  up  time  in  [Q,t]  .  Under  the  conditions  of 
Theorem  3.1, 


yw,  jjpvwdM  «a 
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t~V1[U(t)  -  irQt] 


is  asymptotically 

(3.6) 


(t  -*•  ») 


2 

u 


N(0,1)  where 


.1=1 


2  -1 
Vi 


3.5  Corollary: 

Let  D(t)  be  the  system  down  time  in  [0,t]  so  that  U(t)  +  D(t)  =  t  . 
Under  the  conditions  of  Theorem  3.1 


t“V1[D(t)  -  (1  -  it0)tj 


and 


t"5S[Timi1  +  ViWi“i3]  ^Di^^  ”  iritJ 

are  asymptotically  (t  -*•  «)  N(0,1)  where  D^(t)  is  the  down  time  of  component 

i  in  [0,t]  . 

Applications  of  these  results  to  the  problem  of  determining  maintenance  costs 
are  discussed  in  Section  4. 

Proofs  of  Theorems: 

For  mathematical  convenience,  we  will  work  with  the  random  variable  which 
is  the  number  of  completed  repairs  of  component  i  in  [0,t]  .  We  do  this  because 
we  want  to  assume  that  all  c< — ponents  are  new  at  t  ■  0  .  A  natural  cycle,  then, 
ends  with  a  completed  repair.  Asymptotically,  the  number  of  failures  properly 
normed  will  have  the  same  distribution  as  the  number  of  completed  repairs 
similarly  normed.  For  this  reason  we  use  N^,(t)  to  mean  the  number  of  completed 
repairs  in  [0,t]  for  the  remainder  of  this  section. 

ffiBSTii  nil  ata  tfn  i  r  feMtsaa  aateaaaaaaaa  g&asta  s*aa aa,  - 


YWf 


— 1 — .— ^  - T-  •  p  TT^  ‘(T  »■  *  •  »#V 


Let  S^n  =  £  Xir  and  S^q  =  0  .  The  time  to  the  first  completed  repair  of 

r=l 

component  i  will  then  be 


Vxn> 

Yil  “  Xil  +  Dil  +  sl  DiS 


00 

where  {N.(x);x  0}  is  the  renewal  counting  process  associated  with  jX  l 

J  *  3r  r=l 

as  in  Section  1.  Similarly,  the  time  between  the  r  -  1  and  rth  completed 


repair  will  be 


(3.7) 


r  VSi,r> 

Yir  =  Xir  +  Dlr  +  l  I  D1s 

ir  ir  ir  tfi  s=N.  (S .  J+l  js 

J 


(3.8) 


Z.  =  l  Y. 
in  ir 

r=l 


n  “j  <Si„> 

Sln  +  l  °ir  +  l  l  E  • 

ln  r=l  lr  j*i  r=l  Jr 


3.6  Lemma: 


If  0  <  ui  <  “  and  0  ±  <  00  <1  =  1,2,  , . . ,  k)  ,  then 


“1-  -1  def 

lim  n  Z,  =  p , it..  =  m,  . 
_  in  a.s.  i  0  i 

n-K» 


Proof : 


-1  -1 

By  the  strong  law  n  S,  -*•  p,  and  n  J  D.  -*■  v.  .  Since 

in  a.s.  i  ir  a.s.  i 

r=t 

P,  >  0  ,  S  ■+  <=°  as  n  -*•  «  ,  and  since  p  >  0  ,  N.  (t)  -*■  «  as  t  -»■  »  . 

Ill  &«St  "  ^ 


y  '  a.s. 
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N.(s.  ) 

1  3  l  0 


8j <V  s-1 


js 


Vis 


In 


-][-?]  ats.  bVi 


again  by  the  strong  law  and  (2.3).  The  conclusion  of  the  lemma  follows. 


The  processes  {Z^n;n  JL  U  and  {N^(t);t  >_  0}  are  related  by 


(3.9) 


N±(t)  = 


max  (n  |  Z^n  t)  ,  <_  t 


Zil  >  * 


Since  the  partial  sum  process  {Z^n;n  —  and  the  counting  process  {N^(t);t  0} 

are  essentially  inverses  of  each  other,  it  will  be  sufficient  to  determine  the 
asymptotic  normality  of  the  partial  sum  process. 

From  (3.9)  we  observe  that  for  fixed  v  >  0 


N^ni) 


(3.10) 


I  !Ylr  -  "i1  i  -  VnT)* 

r=l  j_  i  J 


N^(nt) 


"I1"'15  l  [Yir  -  “l1  + 


r=l 


i.N^nxJ+l 


so  that  asymptotically 


N. (nx) 


(3.11) 


i x  7 

n^fnTm”1  -  N^nt)]  -  m^n-"*  £  [Yir  -  m^,]  , 

L  J  r*=l 


where  -  means  asymptotic  equivalent  in  distribution.  On  the  other  hand,  it  is 
well  known  (Renyi  (1957))  that 


n 


-h 


Ni(nT) 

I 

r=l 


[Yir  -  mi]  -  n 


-h 


r-l 


lYir  '  ml! 


(3.12) 
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N.(nT) 

since  -  -*■  ]/e  by  (2.4). 

nt  a.s.  i 


It  will  be  useful  to  expand  Z^n  =  n  ^  £  [y^r  -  m^]  as  follows: 

r=l 


*  i  i  n 

Z,  =  n^S,  -  ny .)  +  n-'5  £  (D,  -  v.) 

in  in  i  v  ir  i 

r=l 

.  W  k  r 

+  »'S  !  I  tt  -vJ+rT5  J  V  N  (S  ) 

j^i  r=l  3  3  j^i  3|_  3  1 


+  n_!l  l  -3-  (S.  -  ny.) 

"j  in  1 


*  -!< 


Let  S.  =  n  5(S.  -n.)  and  c.  =  1  +  J  v.y.  .  Then  we  can  rewrite  Z.  as 

ir.  in  i'  i  . j.  j  ]  in 


Z,  «  c,S,  +  n~  '  l  (D,  -  v.) 

in  i  in  ** ,  ir  i' 

r*! 


,  Nj(sin>  ,  r 

n  l  l  (Djr  -  V  )  +  n_J5  l  v  N  (S  ) 
j^i  r-1  3r  3  iH  3|_3 


["vj1] 


(3.13) 


*  *  v  L  V  J 

Z.  -  c.S,  +  n  ^  )  y  (D.  - 

in  1  in  j=l  r=l  3r 

+  Ji  “"Sh  (n|,i)  '  "V}1]  • 


Note  that  the  summands  in  (3.13)  are  now  independent. 
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3 . 7  Lemma : 


(3.14) 


11m  Var  Z 
n->* 


-  22,  r  2-1,  r  22-3 

in  "  Vi  +  -i  l  Vi  +  »i  I  Wj 


def  2 

=  . 


where  c. 

l 


=1+1 

j^i  J  3 


Proof : 


Use  representation  (3.13)  and  the  well-known  result 


-h  2  —3 

Var  n  [N^  (nt)  -  m/vu]  -*■  T 


We  use  (3.11)  and  (3.12)  to  write 


(3.15) 


n“S[w  '  s!r  -ii  • 

L  J  j  J 


Proof  of  Theorem  3.1: 


From  Theorem  1.2  and  representation  (3.13)  it  Is  obvious  that  the  marginal 


random  variables  N^(t)  are 


/  2  -3\  2 

asymptotically  N(0,w^m^  j  where  w^  is  given  by 
(3.14).  Using  (3.15)  we  see  that 

k  lnwiN 1 

*  *  r  -1  *  ~h  r  L  rJ  J  , 

-  Vi„  -  i  Yj  s,r  -n +  n  X  1 ,  (V "  V 

j^i  J  J  jlmi^  J  j-1  r«l  J  J 


or 


(3. 16)Z 


in 


-1  *  £  -1  * 

'  *0  Sin  _  *  Vjpj  s,r  -n 
j=l  J  jlmy^  J 


_n  +  n_is  l  l 

1  J-1  r=l 


Kv 


®3r  - 


V 


■fc  •‘J-3SV 


£ 


i 


It 
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From  representation  (3.16)  it  is  easy  to  see  that  aribtrary  linear  combinations 
of  the  Z  's  are  asymptotically  normal.  It  follows  that  N^COjl^Ct),  . ..,  N^(t) 


is  asymptotically  multivariate  normal. 

It  remains  to  compute  the  variance-covariance  matrix.  From  Theorem  1.2  and 
I.emma  3.7,  it  follows  that 


2  -3  2 


v .  .  =  v,  ,=  m.  w, 
ii  i  l  i 


as  given  by  (3.3). 

Nj  (t)  ,i^  (t)j  for  i  ^  j  ,  recall  that  by  (3.15) 


i*(m)  -  . 

KJ 


Hence, 


.  d=*  Cov  ^N^nT)  ,N  (nx)J  =  m^m^  Cov 


rz* 


,Z 


*03 


for  0  £  t  £  1  .  Using  (3.16)  we  see  that 


(3.17)  Z 


-1  * 

'  S  _ 


rnr  ]  0 

h. 


nTv 


r  -1  * 

l  V  p  XS 
u  S  S  *“ 


["’Vs1] 


S  =  1 


nTV 


s  J 


+  n"1*  l  '  l 


s=l  r=l 


:d 

sr 


V 


It  fellows  that  for  i  f  j  , 


B^^Sj^a^f^ssa^^^aai^se^^gaa^aaaiaaaaB^a  -  - 


26 


(4.3) 


lin 

t-V® 
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c2^  i  -i 

t  ’  *0  diVi 


} 

A 

A 


by  (2.7). 


4.1  Theorem: 


Under  the  conditions  of  Theorem  3.1, 


(4.4) 


-  «0  ji  dtvi1]  * 


where 


,2  r  r  2  -1  .  2  -3  21.2 

S2  _  J  [,l“l  +  Vi  wiJdi  • 


where  wi  is  given  in  Lemma  3.7. 


Proof : 


■1 


C2(nt)  -  n™Q 


l  ‘‘Vi1] 


k  Ni(nT) 


-  n  **  l  l  di<Dlr  ~  vi)  +  n_ls  l  V  d  (*N  (nr)  -  nxm"1] 
i-1  r-1  1  lr  1  i-1  1  1L  1  1  J 

k  nTmi1  k 

"  n~  4  l  d  (D  -  v  )  -  l  v  d»;Vf  _1-|  . 

i-1  r-1  x  i-1  ijnT°i  J 


Obviously,  the  normed  cost  function  is  asymptotically  normal.  Its  variance  is 
easily  seen  to  be  (when  t  «  1) 


$2  •  j  +  vAiyi 


T  sr’-TiJ  is 


BSTO» r'~r?y^  Tf!^i^7^^9W>wv«)»  ,»?>fy  w\5'  •„  ’, 
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